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Abstract 

The three-dimensional temporal instability of rotating boundary layer flows is investigated by 
computing classical normal modes as well as by evaluating the transient growth of optimal distur- 
bances. The flows examined are the rotating Blasius (RB) and the rotating asymptotic suction 
layers (RAS), with the rotation axis normal to the basic flow plane. In agreement with an inviscid 
criterion, streamwise unstable modes are found in both flow cases for anti-cyclonic rotation: at 
high Reynolds numbers, one obtains the Rossby number unstable range < 1/Ro < 0.57 for RB, 
or < 1/Ro < 1 for RAS. Critical Reynolds and Rossby numbers are also determined in both 
instances. Moreover the dependence of transient growth with respect to wavenumbers, Rossby and 
Reynolds numbers is presented for both cyclonic and anti-cyclonic regimes. In particular, the peak 
transient growth is computed for a wide range of parameter values within the cyclonic regime and 
is shown to be reduced by rotation. A scaling analysis with respect to the Reynolds number is 
performed showing that the standard Re^ scaling is recovered only at very weak rotation. Opti- 
mal disturbances resemble oblique vortices. At weak rotation, they are almost streamwise vortices 
though their structure departs from the classical non-rotating case. Strong rotation imposes two- 
dimensionality and the optimal disturbances vary weakly in the spanwise direction and exhibit 
growth by the Orr mechanism. 

Keywords: Instabilities: Boundary Layer, Boundary Layers: 3D, Geophysical Flows: Rotating Flows 
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I. INTRODUCTION 



The standard development of boundary layers may be transformed by various effects, 
the modification of boundary conditions (riblets, compliant walls) or the presence of bulk 
forces such as centrifugat^ and/or Coriolis forces,— to name just a few. In astrophysical, 
geophysical and technological flows, rotation through the Coriolis effect plays a major role: 
technological examples appear most notably for turbomachinery but also include rotary 
atomization nozzles.— Astrophysical examples include the solar tachocline^ and, on a larger 
scale, accretion disk flow at the disk-star boundary layer.— 

For shear flows {U{x,y),V{x,y),0), it is known that rotation along the 2;— axis can pro- 
duce an instability generating streamwise oriented vortices, also known as longitudinal rolls.—. 
This instability arises due to the imbalance of pressure and Coriolis forces.-'' Rotating channel 
flows have been shown to exhibit this instability, in both theory and experiment,—!^ while 
rotating free shear flows have been predicted to be likewise unstable.— 

However, this modal instability is present only for a certain range of rotations rates and 
depends on the direction of the frame rotation vector (here the 2;-axis) with respect to the 
vorticity vector of the base flow. When the latter quantities are parallel (anti-parallel), the 
rotation is called cyclonic (anti-cyclonic). Cyclonic rotation and strong anti-cyclonic rotation 
are found to stabilize normal modes. On the contrary, weak anti-cyclonic rotation results in 
unstable streamwise vortices. All these properties have been verified for rotating Poiseuille 
flow,'''^'^ rotating planar wakes and rotating mixing layers. In addition, instability thresh- 
olds have been determined in terms of Rossby and Reynolds numbers (the Reynolds number 
Re and Rossby number Ro are defined explicitly in §III(A)). Such an instability mechanism is 
fundamentally inviscid and is reduced by the presence of viscosity. For instance, in rotating 
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Poiseuille flow," a critical point was found at Rec ~ 89 and Rossby number Ro = 1. Even 
though the boundary layer over a flat plate continues to be an essential prototypical flow in 
the study of transition to turbulence, rotating boundary layers have been less extensively 
studied. In one early study/^ it was found that rotating Blasius flow (denoted hereafter 
by RB) is destabilized by weak anti-cyclonic rotation. However stabilization for stronger 
anti-cyclonic rotation was not found since strong rotation rates were not considered. This 
study was limited to a few fixed values of streamwise wavenumbers and no critical point 
was sought. The first purpose of the present work is to investigate the characteristics of 
the unstable normal modes of rotational boundary layer flows by emphasizing the ranges of 
Reynolds and Rossby numbers where the rotational instability is active. 

In non-rotating shear flows (plane Couette, Poiseuille, boundary layer flows), nearly 
streamwise vortices have also been observed during the transition to turbulence (note that 
no fundamental connection exists between these two kind of streamwise vortices). These are 
thought to be the result of a transient growth process that can occur even when linear sta- 
bility theory predicts that all modes are stable. The ability of disturbances in shear flows to 
grow transiently was first proposed by Kelvini^ and explored by Orr;— later, Ellingson and 
Palm,^"^ Landahl^^ and Hultgren and Gustavssoi>i^ proposed mechanism s for this growth. 
Transient growth theory has taken shape more recently see e.g. papers 
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201 reviews. 



23l For steady flows, transient growth can be related to the 
non-orthogonality, with respect to a suitable energy norm, of the eigenvectors of the linear 
evolution operator. This feature permits certain initial conditions to experience relatively 
large amplification G{t) at intermediate times even for asymptotically stable flows. The most 
amplifled transient disturbances for boundary layersi^ typically take the form of streamwise 
vortices. The well known streaks seen in such flows at high levels of free stream turbulence. 



are now thought to be formed as these optimal vortices hft up^^ low velocity fluid from the 
wall region and push down high velocity fluid toward the wallj^ In rotating boundary layers, 
there exists also a wide range of rotation rates for which the instability is not present and 
little attention has been paid to this rotationally stable regime. In particular, the possibility 
that disturbances experience signiflcant transient growth in rotating boundary layers, has 
not been addressed yet (although some preliminary results for rotating channel flow appears 



in Ref. |22l). The second and major purpose of this paper is to study transient growth in 
rotating boundary layer flows, primarily in the regime where no unstable normal modes are 
present. Transient growth factors, maximized in time, are calculated for a broad range of 
wavenumbers and Rossby numbers. Peak transient growth factors, maximized both in time 
and in wavenumber space, are also computed as a function of rotation. Calculations based 
on RB flow are presented alongside equivalent calculations based on the rotating asymptotic 
suction (denoted hereafter by RAS) proflle. The base flows are discussed in §11, then in 
§111 the disturbance equations and calculation methods are presented. Results for modal 
instability and transient growth for both base flows are contained in §IV and the study is 
concluded in §V. 

II. BASE FLOWS 

To examine the effect of rotation on boundary layers, we choose as base flows two well- 
known boundary layer solutions which are rapidly recalled in the following sub-section: the 
Blasius boundary layer and the asymptotic suction proflle. 
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A. Base Flows without rotation 



The paradigmatic Blasius velocity profile U{y)^ is related to the streamfunction / by 
U{y) = U*f'{r]) where U* denotes the free-stream velocity, primes the differentiation with 
respect to the similarity variable tj. More precisely / and t] satisfy 

/'" + ^ff" = with T] = y/A{x) and A(x) = ^/ux/U* . (1) 

The quantity A{x) denotes a boundary layer thickness defined in terms of the free-stream 
velocity U*, the downstream distance x and the kinematic viscosity u. Equation (jT)) is 
easily evaluated using straightforward numerical schemes (Falkner-Skan routines built into 
Matlae^^ were used in this work). Blasius flow includes as well a weak transverse velocity, 



yiy) = y'^ivf'-f), (2) 

which is zero at the wall and approaches the constant displacement velocity, V* = 
as ?/, r/ ^ oo. The displacement thickness 



Sb = (^l-^^dy = A{x) \im[r] - f{l^)] = 1.7208 A{x) 



(3) 



defines the length scale used to build a Reynolds number Re = U*6b/i^- 

The non-rotating asymptotic suction profile (denoted hereafter by AS) describes the 
velocity field obtained when a fluid flows parallel to a surface on which uniform suction V* 
is applied. This boundary layer profile 

U{y) = U*{l-e-y/'^); V{y) = -V* with 6s = ^ (4) 

is a true Navier-Stokes solution and, contrary to Blasius flow, does not vary in the down- 
stream direction. The linea r^^i^^ and nonlinear— stability properties of flow (@]) have been 



long known: AS flow remains stable for much higher Reynolds numbers than Blasius flow 
for a Reynolds number based on the displacement thickness (for AS equal to 6s) and free 
stream velocity, i.e. Re = U*6s/i^- Note that the non-dimensional transverse velocity is 
such that V*/U* = 1/Re. 

B. Base Flows with rotation 

Let us now consider a frame rotating along the z-direction with a constant angular velocity 

fl (see Fig.l). By applying the usual boundary layer scalings^^ to the rotating Navier-Stokes 

equations, the nondimensional equations read in the limit of large Re, 

du ^dv _ ^ 
dx dy 



du du dp d'^u 1 . . 

dx dy dx dy"^ ReRo 



° = -a^-R^" 

where the Rossby number is defined such as Ro = U* /2Q6 where 6 stands for the displace- 
ment thickness {6 = 6b or 6s)- Unless rotation is very strong (i.e. the rotation number 
Rot = (2Ro)-^ » 1), the term V is weak relative to the leading order terms (the 
strongest rotation numbers used in this study are C(l)) and can thus be neglected in the 
streamwise momentum equation, (jH)). A non-rotating boundary layer profile such as the Bla- 
sius profile or AS profile, is then preserved in the presence of rotation^^ although rotation 
calls for an added cross-stream base pressure gradient, dP/dy = —^^U^y) to balance the 
Coriolis force resulting from the streamwise velocity U{y). The rotating asymptotic suction 
(RAS) profile continues to be an exact solution of the rotating Navier-Stokes equations if, 



in addition to the normal base pressure gradient, dP/dy — —^U{y), a streamwise pressure 
gradient dP/dx = ^;^V is included in the base flow to balance the Coriolis term in the 
streamwise momentum equation. 

In the following stability and transient growth calculations, the downstream variation 
of boundary layers is not accounted for. Two major hypothesis follow: (a) the streamwise 
profile is considered at a given streamwise position x — Xq, and (b) the transverse velocity 
component is discarded, i.e. we set V — 0. For all RB computations and most RAS 
calculations (unless otherwise stated), the approximated profile (C/(y),0, 0) will be hence 
used. However, as mentioned above, RAS flow with y 7^ is an exact parallel solution 
of the Navier-Stokes equations; it can thus be used as a test for the validity of the second 
hypothesis. In the Appendix, a comparison between the two possibilities V ^ and = 
is presented for RAS, based on growth rates and neutral curves of normal modes as well as 
transient growth. This study shows that the effect of the transverse velocity on the normal 
modes and transient properties is negligible. 

III. MATHEMATICAL FORMULATION 
A. Linearized disturbance equations 

In a frame rotating along the ^-direction with a constant angular velocity fl (see Fig.l), 
the dimensionless Navier-Stokes Equations read 

— ^ + Utot ■ Vutot + —z X utot = -Vptot + — V^Utot; V • u^t = (8) 
dt Ro Re 

where non-dimensionalization has been performed using free stream velocity U*, displace- 
ment thickness S (= Sb or 63), and pressure scale p{U*Y . Two dimensionless number then 
arise: the Reynolds number Y{b — U*5/v , and the Rossby number Ro = U*/2Q6. 



Following standard methods, the base state velocitj^ iU{y), V, 0) and pressure P{y) are 
perturbed by adding an infinitesimal disturbance {up, Vp, Wp) and Pp. These linear perturba- 
tions can be assumed to satisfy a modal behavior in x,z, viz. 

where a, (3 are assumed to be purely real. By using normal velocity v{y,t) and normal 
vorticity 

C{y,t) = iPu — iaw, (10) 
the linear dynamics is described by one equation governing normal velocity v: 

+ iaU + VD)(D^ - k^)v - ia(D^U)v + ^i(3C - ^(D^ - k'^fv = 0, (11) 
ot Ro Re 

and a second equation governing normal vorticity (: 

(I- + + VD)( - ^iPv + zP{DU)v -^{D'- k^)( = 0, (12) 
ot Ro Re 

where k'^ = D = d/dy. Equation ()11|) (resp. p2|) ) is nothing but the Orr-Sommerfeld 

(resp. Squire) equation modified to account for rotation. Boundary conditions 

v = Dv = C = Q (13) 

should be imposed at the wall (y = 0) and in the far field (|/ — > oo). A dynamical system is 
thus defined for q(y, t) = {v{y,t)X{y,t))'^- 

B. Normal mode approach 

In this work, we study the above problem using two different approaches: the normal 
mode approach (cf asymptotic dynamics) and the non-normal approach (cf transient dy- 
namics). The first approach consists in recasting the problem in §III(A) into an eigenvalue 



problem by assuming q = q(y)e where q = (f , C)^ and uj = ujr+iuJi are complex. Clearly 
a mode having tUj > is unstable. Equations (fTT|) - (fT^ can then be written 



ituMq = Lq 



(14) 



where 



with 



M 



1 



^os T^os 



Cos = {iaU + VD){D^ - k^) - ia{D^U) - ^{D^ - k^f, 

Re 



(15) 



(16) 



CsQ = laU + VD - —{D'' - k^ 
lie 



(17) 



(18) 



It should be emphasized that as a result of rotation, the Orr-Sommerfeld equation is ex- 
plicitly coupled to the Squire equation since both TZqs and IZsq are nonzero while only 
the Squire equation is coupled through TZsq in the absence of rotation. These equations 
together with the boundary conditions 



(19) 



imposed at the wall {y = 0) and in the far field {y oo) define the eigenvalue problem 
to be solved. A Chebyshev collocation algorithm developed in previous work^° is imple- 
mented. First the semi- infinite domain is approximated by a sufficiently large interval [0, L], 
a method known to be effective^S'. In practice, the domain extended to sixteen displacement 
thicknesses 6, L = 16 6. This spatial interval is then mapped to the interval [—1, 1]. Finally 
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eigenfunctions are expanded into Chebyshev polynomials. The number of polynomials 
was generally = 75, although as many as = 125 were used to give improved accuracy 
in some calculations. Collocation points are the standard Gauss-Lobatto points. The eigen- 
value problem (jl4|) and boundary conditions (|19|) are then transformed into a corresponding 
2N X 2N matrix problem. To facilitate calculation, the eigenvalue problem (|T^ is actually 
written as 

ia;q = £q with £ = M^^L . (20) 



C. Transient Growth Calculations 

Transient amplification calculations have also been implemented on the problem in 
pII(A) according to the singular value decomposition (SVD) method given by Reddy and 



Henningson^S. and the algorithms published in Ref. |23l The main steps of the method are 
briefly described below; more details can be found in the references. First, the standard 
kinetic energy norm 

l|q(t)lll= l\\Dv\' + kM' + mdy (21) 

^0 



is used to measure the magnitude of a disturbance q(t) at time t. At time t, one can define 

G{t)= sup ^S = IK"*III, (22) 

where q(0) is an initial disturbance and C is the linear operator defined in (j2nj. This 
quantity G{t) represents the maximum possible energy amplification at time t, optimized 
over all possible initial disturbances. Finally the maximum or optimal growth is defined as 
Gmax = supj>Q G{t). The peak value Gp{Ro, Re) = sup^,^^ Gmax{o:, (3, Ro, Re) can also be 
computed. Note that Gmax is associated with a particular initial disturbance which reaches 
Gmax at a specific time tMAx{c(, P, Ro, Re) such that Gmax = G{tMAx)- Similarly one may 
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define a specific time tp(Ro, Re) related to Gp. 

The quantity G{t) is obtained as follows. Assume that the eigenvalues of (pljl -lfT ^ are 
sorted in order of decreasing imaginary part uOi. One can approximate G{t) by computing 
the maximum possible energy amplification at time t, optimized over all possible initial 
combinations of the K eigenf unctions associated with the first K eigenvalues of ()20|) . In 
this case, it is possible to transform the energy norm of the matrix exponential ()22|) to an 
ordinary 2-norm via the following equation 

||e'^*|||^ ||Fe-'^*F-i|p, (23) 

where A is a x i^' diagonal matrix with the first K eigenvalues uji,uj2i ■■■i^k along the 
diagonal and F is obtained by Cholesky factorization of K x K Hermitian matrix Ai = F^F 
which is calculated using the inner product of the K eigenfunctions of ()20|). the inner product 
used being the same as the one defining the norm (j^D). 

The approximation of (j23p is a result of the finite number K used in the expansion. Fol- 



lowing Ref 



31 



K is chosen such that convergence is reached; in practice, K ^ 60 was found 
to be sufficient. The 2-norm on the RHS of ()23|) can now be evaluated using SVD. This 
procedure gives both Gmax and the K expansion coefficients of the disturbance associated 
with Gmax-^- Note that, since eigenfunctions have been expressed as an expansion in 
Chebyshev polynomials, all calculations are ultimately performed in terms of the Cheby- 
shev expansion coefficients. In particular, the energy norm ()21|) and the weight matrix Ai 
are easily recast in terms of the expansion coefficients using the properties of Chebyshev 



polynomials and their derivatives (for details see Ref. 
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IV. RESULTS 



A. Modal Stability Properties 

It can be shown that within inviscid theory,i^ the stabihty of a rotating shear flow is 
analogous to that of Boussinesq convection in which the total base vorticity {2Q — dU / dy) 
plays the role of temperature gradient. The dimensionless instability criterion then states 
that total vorticity must be negative somewhere in the flow^ i.e. 

Ro ^ > 1 . (24) 

dy 

The relative vorticity of the base flow —dU/dy which is directed along the z direction, is 
always negative for both RB and RAS flows. Positive Rossby Ro > 0, thus corresponds to 
oppositely directed rotation and flow vorticity, that is why it is referred to, in the geophysical 
literature, as anti-cyclonic. This is the region we will consider here since negative Rossby 
Ro < correspond to cyclonic case where this instability disappears. 

Neutral stability curves of RB flow at Re = 600 and several Rossby numbers are shown in 
Fig. 2a. The region identified as "TS" corresponds to the (rotationally modified) ToUmien- 
Schlichting modes of the non-rotating Blasius flow. The TS region is only slightly modified 
by rotation, becoming discernibly larger only when Ro is below 300. The lower curves in Fig. 
2a, identified by their Ro values, delineate the region of rotational instability. The largest 
growth rates are found for streamwise uniform modes (a = 0) having a non-zero spanwise 
wavenumber (3 (Fig. 2b), that depends on the value of Ro, in agreement with previous 
results for other rotating shear flows. The extent of the rotationally unstable region 
grows rapidly with rotation, and is quite large at Ro = 300 which is the strongest rotation 
depicted in figure Fig. 2a. Ultimately, this instability is suppressed for very strong rotation 
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as examined below. Note that the rotational instability takes on much larger growth rates 
than the TS instability, as illustrated by the curves of constant growth rate shown in Fig. 
2b for the case Re = 600, Ro = 300. More precisely the largest TS growth rate is found in 
this case to be sup^ p{uJi)Ts = 0.006 while the maximum value in Fig.2b is sup^^^cUj = 0.145 

As in RB flow, an anti-cyclonic rotation for RAS leads, if not too strong, to a rotationally- 
induced instability which also reaches its maximum at a = 0, while a cyclonic rotation 
(Ro < 0) is stabilizing. These rotationally unstable regions in (a,/?)-plane are presented 
in Fig. 3a for Re = 600 and several Rossby numbers Ro. Note that growth rates of the 
rotational instability for RAS are of comparable magnitude to those found for RB flow (see 
Fig. 3b for the case Ro = 300 where sup^pUi = 0.144 is found). On the contrary, no TS 
neutral curves appear for RAS (Fig. 3a) for Re = 600. Recall that the AS profile is known 
to be much more stable with respect to TS instability compared with non-rotating Blasius 
flow. This effect is very apparent in the large critical Reynolds number of the AS profile 
(Rec = 54370) RAS flow thus lacks TS instability at Re = 600. 

Re- stabilization - 

In RB flow, the condition for instability (|24|) corresponds to 1/Ro < 0.5714^^ while in 
RAS flow the condition becomes 1/Ro < 1. The re-stabilization of rotating boundary layer 
flows by strong anti-cyclonic rotation is confirmed by numerical computations. The curves 
in Fig. 4 represent the largest growth rates max uJi of this rotational instability found in the 
(a,/3)-plane as a function of 1/Ro (this maximum was always found to occur for a = 0). 
At the largest Reynolds number depicted in Fig. 4a, Re = 6000, re-stabilization of RB 
flow occurs at 1/Ro ~ 0.57 while for RAS, re- stabilization occurs at 1/Ro ~ 0.98 when 
Re = 60000, in agreement with the predicted values. Large positive rotations have also been 
found to re-stabilize channel flows,^^^ at Ro values consistent with the vorticity criterion 
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(El- 

Since the rotational instability is fundamentally inviscid, the minimum Reynolds number 
below which the rotational instability is suppressed is quite low. The minimum Reynolds 
numbers below which no instability is found are Rec ~ 21 at Roc ~ 10.9 for RB (see Fig. 
4a) flow and Rec ~ 25.9 at Roc ~ 23.5 for RAS flow (see Fig. 4b). 

B. Transient Amplification 

The effect of rotation on transient growth for the RB and RAS flows is quantified through 
the analysis of amplification gain GMAx{(y-i Ro? Re) as a function of streamwise wavenum- 
ber a, spanwise wavenumber /3, Reynolds number and Rossby number. Note that, in calcu- 
lating Gmax and Gp values, a fixed time interval [0, T] is used during which the maximum 
growth is allowed to occur; the value of T required to capture the maximum varies with Re 
and is adjusted accordingly in the calculations. 

Non-rotating boundary layer flows are known to display the following transient growth 
features: (i) the most amplified disturbances tend to be three-dimensional with (3^0 and 
a = 0; (ii) the peak amplification for large Reynolds number is proportional to the square 
of the Reynolds number i.e. 

Gmax oc — , [Zb) 

where it is recalled that fc^ = + /J^ and T is a function of k and aRe only; (iii) the time 
to achieve peak amplification increases linearly with the Reynolds number i.e. tMAx oc Re. 
For rotating boundary layers, scaling ipHjl or its straightforward consequence Gp oc Re^ is 
not preserved as shown below. 

For non- rotating boundary layers {U{y), 0, 0), the Reynolds scalings rely on the following 
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statements (see also Ref|23|): upon performing the rescalings t = t/Re, ( = C/(/3Re), v = v, 



equations (fTT|l - (|T^ become dependent only on k'^ and aRe and the energy turns into 

l|q(t)|||= [\\Dv\' + k'\v\')dy + iPRef !\\C?)dy . (26) 
Jo Jo 

When the Rossby number is present, the energy expression is still valid but, in addition to 

= a'^ + and aRe, the parameters Ro and explicitly appear in equations (fTT|) - (fT^ 

d 1 

+ iaReU)(D^ - P)v - iaRe(D^U)v + —i(pRe?( - (D^ - k^fv = 0, (27) 
dt Ro 

+ iaReU)C - ^^iv + i{DU)v - {D^ - P)( = 0. (28) 
dt Ro 

The standard arguments are thus no longer valid when the Rossby number is present, im- 
plying that the classical scalings will not be replicated. The numerical computations which 
confirm this idea, are presented below with the observed scalings. 

Rotating Blasius boundary layer - 

In Fig. 5, level curves of Gmax in the (a, /3)-plane are presented for RB flow in three cases: 
weak cyclonic rotation Ro < (Fig. 5a), no rotation (Fig. 5b), and weak anticyclonic rotation 
Ro > (Fig.5c), the Reynolds number being always fixed at Re = 600 which is just above 
the critical value Rec ~ 520 for unstable ToUmien-Schlichting waves. The TS instability 
region is apparent along the a-axis centered at a ~ 0.3 and has only minor influence on the 
transient growth (see below). 

In the absence of rotation, the peak Gmax value was found to occur at a = and 
j3 = 0.651 (Fig. 5b), and to take the value Gp = 545.8 ( Gp = 545.9 with N = 75 polynomials 
and Gp = 545.8 with = 125). This value is in good agreement with other studies . 
From previous computations by Butler and Farrell^^ {Gp = 1514 at Re = 1000) or Corbett 
and BottaroM {Gp = 1515 at Re = 1000) the scaling law Gp = 1.515 x lO^^Re^ should 
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be the correct one for non-rotating Blasius flow. Using this latter formula, it is found that 
Gp = 545.4 at Re = 600. 

When the rotation is anti-cyclonic (Fig.5c) i.e. for Ro > 0, rotational instability of §IV(A) 
appears along the /5-axis and the level curves of Gmax are distorted in the neighborhood of 
the highhghted region of unstable normal modes. Moreover, the Gmax values adjacent to 
the neutral curve are enhanced by the nearby modal instability, and the peak amplification 
occurs along the rotational neutral curve at /3 ~ 1. In such a case, transient growth must 
compete with the strong modal instability. 

With cyclonic rotation (Fig.5a) i.e. for Ro < 0, Gmax values are reduced in magnitude 
everywhere in the (a, /3) -plane and the rotational instability is absent. In addition, Gp is 
now found at non-zero streamwise wavenumber (for Ro = —1800 and Re = 600, a ~ 0.12), 
while the spanwise wavenumber slightly shifts from the value of the non-rotating case (for 
Ro = —1800 and Re = 600, (3 ~ 0.75). An optimal disturbance (hereafter denoted by OD) 
with weak streamwise dependence {a ^ 0) is often referred to as oblique, and is more typical 
of plane Couette flow than of boundary layer flow. Here, the obliquity is a consequence of 
the Coriolis effects. The Ro dependence of the peak growth Gp is illustrated in Fig. 6. Each 
point in this flgure was obtained by flrst performing a calculation of Gmax throughout the 
{a, /?)-plane (like the one of Fig.5a), then flnding the peak value, Gp, and the wavenumbers, 
ap and (3p at which the peak occurs. Two Reynolds numbers are shown: Re = 600, marked 
by open circles, and Re — 1800, marked by dots. 

The inset of Fig. 6 shows the location in the a — (5 wavenumber plane corresponding to 
the plotted Gp values. When rotation is very weak, Gp is found at ap = 0, Pp = 0.651; as 
rotation increases, the location of Gp loops through the plane, eventually settling at /3p = 
for strong rotation, with a streamwise wavenumber Q!p(Re) dependent on the Reynolds 
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number. The value Gp = 16.70 is found at ap = 0.467 for Re = 600 while Gp = 51.02 is 
found at ap = 0.364 when Re = 1800. These values are also recovered for strong values 
of anti-cyclonic rotation. The wavenumber ap and Gp values at strong rotation agree with 
two-dimensional non-rotating transient growth calculations. Butler & FarrelU^ give Gp = 28 
at ap = 0.42 for the best optimal 2D {j3p = 0) disturbance for non-rotating Blasius flow at 
Re = 1000. At the same Reynolds number, we find Gp = 28.4 for Ro = —0.5 at ap = 0.42. 
These findings are consistent with the idea that strong rotation imposes a two-dimensional 
Taylor-Proudman regime.^^*^ 

The ratio of the two Gp curves in Fig. 6 is a measure of the Re dependence of Gp over 
a range of Ro. From the figure, this dependence is clearly inconsistent with a Re^ scal- 
ing throughout the entire Ro range. This Re dependence is explicitly shown in Fig. 7a and 
7c. For weak rotation (Ro = —10^), Gp/Re^ possesses a weak Re dependence (Fig. 7a) 
which strengthens for slightly stronger rotations (Ro = —10^). In the strong rotation 
regime (Fig. 7c), no simple Re dependence is found for Gp in RB flow: although Gp seems 
linear in Re, a careful analysis reveals that this is not strictly correct. In addition, figure 
Fig. 7c confirms that, for transient growth properties, the fiow becomes increasingly two- 
dimensional in character as rotation increases. Note that, in a similar manner, the time tp is 
not linear in Re for weak rotation (Ro = —10^) (Fig. 7b) while we obtain scaling tp = 0.78Re 
with no rotation. In the strong rotation regime (Fig.7d), tp possesses no obvious scaling 
though one would guess by eye a V^Re law, which is not confirmed by a rigorous fit. 

Rotating asymptotic suction flow - 

The qualitative features found for RB fiow are also recovered for RAS fiow; the following 
presentation can therefore be relatively terse. Figure 8 is the equivalent of Fig. 5, but 
for RAS fiow at Re = 600. AS fiow is well below its critical value Rc = 54370 for TS 
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instability, so no TS region is present. The most amplified disturbances in the absence of 
rotation (Fig. 8b) are again streamwise uniform {a = 0) while for weak cyclonic rotation 
(Fig. 8a) they begin to show weak streamwise variability {a ~ 0.07 for the case of Fig. 8a). 
Figure 9 is the RAS flow equivalent of Fig. 6; here, again, both Re = 600, marked by 
open circles, and Re = 1800, marked by dots, are calculated. These are the only known 
calculations of transient growth in RAS flow. For AS flow (i.e. no- rotating case), the values 



of Gp and (3p have appeared in recent works: Ref. |2a provides a scaling law of the form 
Gp = 0.99-Re^ which gives e.g. Gp = 356.4 at Re = 600. In the limit of very weak 
rotation, we find Gp(Re = 600) = 358.12 (or, when the transverse velocity V is included, 



36|). As in RB, scalings Gp oc Re^ 



G'p(Re = 600) = 356.74, in good agreement with Ref 
and tp oc Re are not recovered, except for very weak rotation (Fig. 10a and 10b). Again 
as in RB, a careful analysis does not reveal a clear scaling for Gp and tp in the strong 
rotation regime (Fig. 10c and lOd) though one would guess by eye a scaling in Re and v^Re 
respectively. 

The inset of Fig. 9 shows the location in the a — (3 wavenumber plane corresponding to 
the plotted Gp values. As in RB, the optimal disturbance in RAS flow makes an excursion 
in wavenumbers, starting at ap = 0, /3p = 0.499 {jSp = 0.530 when V is included) for very 
weak rotation, looping through the plane and settling on /? = and a Re-dependent ap(Re) 
when rotation is strong. For instance, the value Gp = 5.736 is found at ap = 0.670 for 
Re = 600 while Gp = 13.12 is found at ap = 0.617 for Re = 1800. When V is included, 
these values shift to Gp = 5.661 at ap = 0.675 for Re = 600 and Gp = 13.01 at ap = 0.621 
for Re = 1800. These same values are also recovered for strong values of anti-cyclonic 
rotation. 
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C. Optimal Disturbances 

The most amplified disturbances in cyclonic RB and RAS flows possess a non-null stream- 
wise wavenumber: they thus correspond to oblique vortices. The flow fleld of these dis- 
turbances can be calculated from its expansion coefficients, already obtained in the SVD 
solution of Gmax- An example of the optimal disturbance (OD) in a weakly rotating RB 
flow, at Re = 600, Ro = — 1 x 10^, is depicted in Fig. 11 and can be compared to the OD for 
the classical Blasius flow at the same Reynolds. In flgures Fig. 11 and 12, the velocity fleld 
vy + wz and the streamwise velocity component u are shown at a given streamwise location 
{x — 0). Due to the similarity of the optimal disturbances for RB and RAS flows, RAS 
examples are not shown. 

While the OD in RB at very weak rotation Ro = — 1 x 10^ is indistinguishable from the OD 
in non- rotating Blasius flow at both t — and t — tp, this is not so for larger but still weak 
rotations. Indeed the OD in RB at Ro = — 1 x 10^ (Fig. 11 a) -which resembles with apparent 
counter-rotating streamwise rolls, the OD in non-rotating Blasius flow a,t t — (Fig. 11c)- is 
such that the extrema of streamwise velocity u (Fig. lib) are not precisely aligned with the 
up- and down-flows (Fig.lla) between the vortices, as found in non-rotating Blasius flow 
(Figs. 11c, d). This offset is still found in the OD at its greatest ampliflcation at t = tp. 
Except for a shift of phase, little change is seen in the u fleld. This explains why the Gp 
value is not too different from the one found in the non-rotating limit. The kinetic energy 
of the OD is initially mainly in the v ot w components and not in u component. At t — tp 
it is contained mainly in u fleld. 

The above features are enhanced at stronger rotations as depicted in Fig. 12, where Ro = 
—1800: even the OD at t — being profoundly modifled. Finally, at very strong rotation 
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(e.g. Ro = — 1 in Fig. 13), two-dimensionality is imposed in the x — y plane and the OD 
at t = (Fig. 13) clearly exhibits the characteristic tilting against the shear of the 2D Orr 
mechanism. At t — tp, the shear has evolved the initial disturbance into vortices counter- 
rotating along the 2;-axis (Fig. 13). 

V. CONCLUSIONS 

This work has expanded on previous studies of boundary layer flows to account for frame 
rotation. In particular, the temporal evolution of three-dimensional disturbances has been 
studied for RB and RAS boundary layer flows in the presence of spanwise rotation. Normal 
mode and transient growth disturbances properties were calculated and compared for the 
two flows. Rotational modal instability, similar to that known for rotating channel and 
free shear flows, was found for weak anti-cyclonic rotation in RB and RAS; neutral curves, 
growth rates and the critical point were calculated for both flows. Rotating boundary layer 
flows were also found to exhibit transient growth. In the modally stable cyclonic regime 
the transient growth properties strongly depend on the rotation rate. In particular: (i) 
peak growth factors are reduced by rotation; (ii) even at weak rotation, the classical Re^ 
dependence is not recovered; (iii) optimal disturbances take the form of oblique vortices 
which transform continuously from streamwise uniformity {a — 0) when rotation is very 
weak to spanwise uniformity (/? = 0) when rotation is strong, this last feature being a 
consequence of a Taylor-Proudman regime. 

It is anticipated that the results presented here may be useful in dealing with geophys- 
ical and astrophysical flows, where rotation is important and curvature is often negligible. 
Most noteworthy is the example of Keplerian flow in an accretion disk, where no purely 
hydrodynamic inviscid linear eigenvalue instability has been identifled. A study of transient 
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growth in two-dimensional accretion disks^S has recently been performed and further inquiry 
is already under way. The problems analyzed in this work constitute an interesting element 
of disk flow studies. 

Finally let us mention that the longitudinal roll instability in rotating shear flow exhibits 
weak downstream developmenlj^. Because the development is weak and because a great deal 
of closely related previous work2*2*i24iii2& is based on temporal theory, this study adopted 
a temporal framework. However recent spatial analyses of ordinary boundary layer tran- 
sient growtb22i^Mi are shown to be more realistic with respect to observations.— In future 
investigation, this aspect should be envisaged to complement the present temporal theory. 
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APPENDIX: The transverse velocity in the presence of rotation. 

Transverse velocity V is normally neglected in the classical temporal stability analysis of 
Blasius flow, being of the same order as other neglected terms. The effect of this assumption 



on the instability normal modes is we, 
(see standard references, e.g. Ref 



1 documented for non-rotating AS and Blasius flows 
if). In this appendix, we examine the effect of this 
assumption on two aspects which have not been previously studied: (a) the effect of V on 
transient growth properties with or without rotation; (b) the effect of V on the rotational 
normal mode instability properties. We consider these two effects on the AS and RAS flows 
since they are true Navier-Stokes solution when V is present which authorizes a quantitative 
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comparison between the simplified analysis {V = 0) and the complete one. 

An overview of the effect of V on the transient growth properties of the AS flow (i.e. 
non-rotating) can be found in Table I. Briefly, the neglect of V has little effect on the peak 
growth factors Gp of AS flow, less than one percent. The wavenumber of peak growth, (3p, 
is shifted to larger values when V is included in AS calculations, but the shift is relatively 
small. The two cases are illustrated in Fig. 14 where level surfaces of Gp{a, (3) are shown for 
the AS calculated with and without the V term. As regard tp values, we obtain tp = 1.145Re 



when V is included which agrees with RefJ36l and we find tp = 1.331Re when V is neglected. 

The effect of V on the rotational instability neutral curves and growth rates is small, ex- 
cept near the critical point (Rec, Roc). These properties are illustrated for RAS in Fig. 15a, 
showing the shifted neutral curves, and Fig. 15b, showing the modified growth rates. When 
V is included, the location of the critical point shifts from Rec = 25.9, Roc = 23.5 to 
Rec = 33.7, Roc = 14.6, as seen in Fig. 16 (compare to Fig. 4(b)). 

Under rotation, the effect of V on the transient growth properties of the RAS remains 
negligible. To show this we display in Fig. 17 the differences in peak growth factors Gp 
calculated with and without the V term for Re = 600, Ro = 1800. 
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TABLE I: Peak transient growth quantities for non-rotating AS flows computed with and without 
the transverse velocity, V term. Values for non-rotating Blasius flow with V = are given for the 
sake of completeness. 



Blasius Blasius AS AS 

Re = 600 Re = 1800 Re = 600 Re = 1800 

Gp{V = 0) 545.93 4910.5 358.12 3220.2 

Gp{V^O) — — 356.74 3207.5 

/3p(V=0) 0.651 0.651 0.499 0.499 

Pp{V^O) — — 0.529 0.530 

Table I, Yecko, Phys. Fluids 
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List of Figure Captions 

FIG. 1: Sketch of the flow configuration. 

FIG. 2: (a) Neutral curves of RB flow at Reynolds number Re = 600 and various Rossby 

numbers Ro; a and /? are the streamwise and spanwise wavenumber, respectively, (b) Curves 
of constant growth rate loi for RB flow at Re = 600 and Ro = 300. 

FIG. 3: (a) Neutral curves of RAS flow at Reynolds number Re = 600 and various Ro; (b) 
Curves of constant growth rate loi for RAS flow at Re = 600 and Ro = 300. 

FIG. 4: (a) Largest growth rates max (Vi for anti-cyclonic RB flow at Reynolds numbers 
Re = 25,60,600,6000 (curves bottom to top); Inset: projection of the neutral surface 
uJi{P, Re, Ro) — onto the Re — Ro plane, showing the critical point Rec — 21.005, Roc — 
10.9; (b) Same as (a) but for RAS at Re = 40,60,600,6000,60000 (curves bottom to top); 
Inset: critical point Rec — 25.9, Roc — 23.5. 

FIG. 5: Level curves of Gmax in the (a, /3)-plane for RB flow at Re = 600 and for weak 
cyclonic (Ro = —1800) rotation (a), no rotation (b), and weak anti-cyclonic (Ro = 1800) 
rotation (c). Regions of modal instabihty are shaded. 

FIG. 6: Peak amplification factors Gp of RB fiow at Re = 600 and Re = 1800 as a function 
of cyclonic Ro < 0; Inset: corresponding location in the wavenumber a — (5 plane of the 
optimal disturbance associated with Gp values in the figure; weak rotation cases lie along 
the 13 axis while strong rotation cases he along a. 

FIG. 7: (a) Gp/Re^ of RB flow as a function of Re for two weak rotation cases Ro = —10^, 
Ro = —10^ and no- rotation; (b) tp/Re of RB flow as a function of Re for the same cases as 
(a); (c) Gp of RB flow as a function of Re for strong rotation, Ro = —1, and two-dimensional 
disturbances; (d) tp of RB flow as a function of Re for the same cases as (c). 
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FIG. 8: Level curves of Gmax in the (a, /3)-plane for RAS flow at Re = 600 for weak cyclonic 
(Ro = —1800) rotation (a), no rotation (b), and weak anti-cyclonic (Ro = 1800) rotation 
(c). Regions of modal instability are shaded. 

FIG. 9: Peak amplification factors Gp of RAS flow at Re = 600 and Re = 1800 as a 
function of cyclonic Ro < 0; Inset: corresponding location in the a — (5 plane of the optimal 
disturbance associated with Gp values in the figure; weak rotation cases lie along the (5 axis 
while strong rotation cases lie along the a axis. 

FIG. 10: (a) Gp/Re^ of RAS flow as a function of Re for two weak rotation cases Ro = —10^, 
Ro — —10^ and no-rotation; (b) tp/Re of RAS flow as a function of Re for the same cases 
as (a); (c) Gp oi RAS flow as a function of Re for strong rotation, Ro = —1, and two- 
dimensional disturbances; (d) of RAS flow as a function of Re for the same cases as 
(c). 

FIG. 11: (a) Velocity vectors vy + wz of the most amplified disturbance in cyclonic RB flow 
at Re = 600 and Ro = —10^. The disturbance is shown at t = (left) and at the time 
of maximum growth (right), (b) Streamwise velocity u for OD of (a), (c) Velocity vectors 
vy + wz for non- rotating Blasius flow at Re = 600. The disturbance is shown at t — (left) 
and at the time of maximum growth (right), (d) Streamwise velocity u for OD of (c). 

FIG. 12: (a) Velocity vectors vy + wz of the most amplified disturbance in cyclonic RB flow 
at Re = 600 and Ro = —1800. The disturbance is shown at t = (left) and at the time of 
maximum growth (right), (b) Streamwise velocity u for OD of (a). 

FIG. 13: Velocity vectors ux + vy in the x — y plane of the most amplified disturbance in 
cyclonic RB flow at Re = 600 and Ro = —1. The disturbance is shown ai t — (left) and 
at the time of maximum growth (right). 
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FIG. 14: Comparison of the level curves of Gmax in the {a, /3)-plane for AS flow calculated 
with V included (solid curves) and V neglected (dashed). The peak growth is found at 
ap = and Pp = 0.529 (V included) or /3p = 0.499 (V neglected). 

FIG. 15: (a) Neutral curves of RAS flow at Reynolds number Re = 600 and three Ro; (b) 
Curves of constant growth rate LOi for RAS flow at Re = 600 and Ro = 300. In the figure, 
the dashed curves correspond to calculations in which V was neglected, reproducing Fig.3; 

the solid curves were computed with V included. 

FIG. 16: Largest growth rates max coi for anti-cyclonic RAS flow with V included at Reynolds 
numbers Re — 40,60,600,6000,60000 (curves, bottom to top); Inset: projection of the 
neutral surface u!i{/3, Re, Ro) — onto the Re — Ro plane, showing the critical point at 
Re = 33.77, Ro = 14.5. To be compared to Fig.4(b). 

FIG. 17: Comparison of maximum transient growth level curves for RAS flow at Re = 
600, Ro = 1800, calculated with V included (sohd curves) and V neglected (dashed curves). 
The peak growth factor is Gp = 130.79 {V included) or Gp = 129.86 {V neglected). 
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